REVIEW 4 major objections 4 minor 49 references
No single solar latitude follows the sunspot cycle exactly.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
Supergranule lane widths and Ca K intensities correlate most strongly with the sunspot cycle at different latitudes and with different time lags, so no single latitude follows the cycle exactly.
T0 review reviewed 2026-08-04 challenge →
load-bearing objection A useful 100-year latitude–lag map, but the headline intensity lag of 1.25–1.5 yr is likely an interpolation artifact from yearly sampled data. the 4 major comments →
Which Solar Latitude Follows the Sunspot Cycle Exactly?
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The central discovery is that the supergranular network responds to the solar cycle differently in its geometry and its brightness. Lane widths—a proxy for magnetic flux concentrated at cell boundaries—follow the sunspot cycle almost perfectly near ±20° latitude, with no measurable phase delay. Mean network intensity, however, correlates most strongly near ±13–14° and peaks 1.25–1.5 years after sunspot maximum. Thus no unique latitude can be said to follow the sunspot cycle exactly; the correlation surface has broad, symmetric peaks whose latitude and lag depend on the quantity measured. The paper also finds significant, cycle-correlated lane-width variations across 55°S to 55°N, whereas int
What carries the argument
The key quantities are the supergranular network lane width and its mean Ca II K intensity, extracted from small image windows of century-long Ca II K spectroheliograms after equal-contrast calibration, limb-darkening correction, and rejection of the top and bottom 5% of window intensities. Lane width is measured as the width of the autocorrelation function of each window. The argument is carried by latitude-by-latitude cross-correlation functions between these yearly averaged time series and the sunspot number, with the lags refined by interpolating the correlation curves to 0.01-year resolution. This machinery isolates which latitudes lock onto the cycle and at what delay.
Load-bearing premise
The load-bearing premise is that the equal-contrast calibration of the century-long Ca II K series—especially the trial-and-error rejection of the top and bottom 5% of window intensities—preserves true latitude-dependent changes in lane width and intensity; if seeing degradation or the threshold choice biases some latitudes, the correlation peaks and lags could shift.
What would settle it
Compute the same latitude-resolved cross-correlation on raw, uncalibrated Ca II K contrast data and on a modern space-based network-brightness dataset. If the lane-width peak moves away from 18–20° by more than the stated ±2°, or if the intensity peak no longer lags by 1.25–1.5 years, the central claim would be falsified. A simpler check: split the 1907–1990 record into early and late halves; if the latitude and lag of maximum correlation are not stable across halves, the result is likely an artifact of calibration drift.
If this is right
- Lane width at roughly ±20° can serve as a phase-zero tracer of the sunspot cycle, useful for cycle timing and amplitude prediction.
- Intensity lags the cycle by 1.25–1.5 years at ±13–14°, meaning brightness-based quiet-Sun indices will peak after sunspot maximum.
- The different latitude and lag structure implies flux transport redistributes network magnetic field before it is seen in brightness, affecting surface flux-transport models.
- Quiet-Sun UV irradiance reconstructions should use latitude-dependent network intensity with a post-maximum lag rather than sunspot number directly.
- The contrast between the broad 55°N–S significant correlation range for lane width and the narrow 5–35° bands for intensity cautions against using a single latitude to represent the quiet Sun.
Where Pith is reading between the lines
- A natural extension would be to track the same lane-width/intensity lag pattern in modern space-based magnetograms; if the lag reversal near ±25° and ±40° persists, it would strengthen the flux-transport interpretation and yield a direct transport-speed estimate.
- The paper's lane-width/cycle relationship may help resolve contradictory reports on supergranular size versus cycle: if size correlates with lane width only at certain latitudes, comparisons at mismatched latitudes would explain the disagreement.
- The reported difference between zero lag for lane width and 1.25–1.5-year lag for intensity could be tested with other chromospheric or EUV network brightness proxies, where the lag should shift systematically with formation temperature.
- Because the analysis stops before 1990 due to seeing degradation, applying the same method to modern uninterrupted space-based images would show whether the correlation surface is stable across solar cycles 23–25.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes ~34,000 Kodaikanal Ca II K spectroheliograms (1907–1990) to measure yearly averaged supergranular lane widths and intensities as functions of latitude (60°N–60°S). It cross-correlates these quantities with the international sunspot number at each latitude and reports that the maximum correlation for lane widths occurs near 18–20° latitude with zero lag, while intensities peak near 13–14° latitude with a lag of 1.25–1.5 yr after solar maximum. The paper concludes that no single latitude follows the sunspot cycle exactly for all quantities and discusses implications for flux transport and quiet-Sun UV irradiance. The main result—that correlation with the solar cycle is strongest in low-to-mid latitudes—is visible in the figures and is broadly plausible, but the paper's most specific quantitative claims (the sub-annual intensity lag and the quoted latitude uncertainties) are not supported by the statistical analysis as presented.
Significance. If the central claims can be established, the paper would provide a valuable long-baseline, latitude-resolved characterization of how the supergranular network responds to the solar cycle. Its strengths include the use of a nearly century-long homogeneous archive, an external sunspot-number series rather than a self-referential proxy, and a clear presentation of the latitude-lag structure. The finding that lane width and intensity peak at different latitudes and lags, if robust, would be relevant to flux-transport models and to understanding the quiet-Sun UV irradiance cycle. However, the paper currently overstates the precision of its headline phase lag: the CCF is computed from yearly averages and then interpolated to 0.01-yr intervals, which cannot create the missing sub-annual information. In addition, the significance threshold does not account for autocorrelation or multiple testing across 120 latitudes. These issues do not invalidate the broad correlation pattern, but they do affect the quantitative conclusions and need to be addressed before the paper can be accepted.
major comments (4)
- [Section 3, Figure 2 and interpolation paragraph] The headline intensity lag of 1.25–1.5 yr is not established by the analysis described. The cross-correlation is computed from yearly averaged quantities, so it is defined only at integer-year lags. Interpolating the ±5-yr window with 0.01-yr spacing cannot create genuine sub-annual phase information; the reported lag lies between lags 1 and 2, precisely where interpolation of a broad, noisy CCF is most uncertain. The authors should either (a) compute the CCF using sub-annually binned sunspot and network data, (b) fit a parametric model to the CCF with an uncertainty estimate, or (c) provide a bootstrap/permutation distribution of the peak lag. Without such a test, the contrast between the zero-lag lane-width result and the 1.25–1.5-yr intensity lag remains an interpolation artifact.
- [Section 3, Figure 3 top and bottom panels] The significance threshold of 0.29 (claimed 99% confidence) does not account for the strong autocorrelation in the annual sunspot number and in the yearly averaged lane-width/intensity series. The effective number of independent samples is far below the number of years, and scanning 120 latitudes introduces a multiple-testing problem that can easily produce spurious peaks above 0.29. The claim that correlations are 'highly significant' over a broad latitude range therefore needs support from a test that preserves temporal autocorrelation (e.g., block bootstrap, phase scrambling, or ARMA-based effective degrees of freedom) and controls the false-discovery rate across latitudes.
- [Section 2, intensity rejection threshold] The 5% top/bottom intensity rejection threshold is chosen by trial and error, and the manuscript asserts that 'varying the threshold by a few percentage does not significantly affect the correlation' without showing any quantitative test. Because this threshold directly determines which windows contribute to the latitude-dependent lane-width and intensity time series, a few-percent change could shift the correlation peaks and lags that are central to the conclusions. Please provide a sensitivity analysis (e.g., thresholds of 3%, 5%, 7%, 10%) and report the resulting peak latitudes and lags, or at least show that the headline values are stable.
- [Section 3, Figures 2 and 3] Correlation coefficients are reported without uncertainties, and the quoted latitude uncertainties of ±2° have no stated derivation. The text does not explain how the peak latitude or its uncertainty was obtained from the smoothed curves, nor does it give confidence intervals on the correlation values. Since the paper compares the latitude peaks of lane width and intensity (18°N/20°S vs. 13°N/14°S) and interprets their difference, a bootstrap or Monte Carlo procedure that yields uncertainties on the peak locations and lags is necessary to support those comparisons.
minor comments (4)
- [Title and abstract] There are typographical errors: 'F ollows' in the title, 'intensitiy' in the abstract, and 'coefficients' in several places. Please correct these.
- [Section 2, Data & Analysis] The description of the 'equal-contrast technique' is brief; a reader unfamiliar with Raju (2020) and Singh et al. (2021) will not know how the FWHM contrast target (0.10–0.11) is converted into an intensity scale or whether this calibration is stable across the different photographic emulsions over 100 years. A sentence summarizing the main calibration steps would improve reproducibility.
- [Section 3, lag sign convention] The paper reports lags for lane width and intensity but never explicitly defines the sign convention (e.g., positive lag means quantity follows the sunspot cycle). Adding a sentence explaining the convention would prevent misinterpretation of the values in the bottom panel of Figure 3.
- [References] The reference list contains an incomplete author name for Roudier et al. (2014): 'Roudier, T., Švanda, M., Rieutord, ., et al.' The missing author initial/name should be corrected.
Circularity Check
No significant circularity: lane widths/intensities are cross-correlated with the external SILSO sunspot series, so the central result is anchored outside the paper's own fitted values. Self-citations supply only the measurement pipeline and do not force the found peak latitudes or lags. The 1.25-1.5 yr intensity lag rests on interpolating a yearly-sampled CCF, a resolution limitation, not a circu
full rationale
The derivation chain is: (1) lane widths (autocorrelation FWHM) and mean intensities are measured from Kodaikanal Ca II K windows; (2) yearly averages are cross-correlated with external World Data Center SILSO sunspot numbers; (3) latitudes and lags of the CCF maxima are read off. None of these steps defines the target result in terms of the inputs. No parameter is fitted and then renamed a prediction: the equal-contrast normalization, the 5% window rejection, and the autocorrelation width are fixed preprocessing choices from published work (Singh et al. 2021; Raju 2020; Raju et al. 2023). The paper asserts threshold robustness without a quantitative sensitivity test, which is missing support, not a circular step. Self-citations (Raju 2016/2018/2020/2023) supply the measurement methodology and prior evidence of cycle dependence; they are published, parameter-free methods whose assumptions do not include the present peak latitudes or lags, so they constitute independent support and do not force the result. Two flagged limitations are correctness risks, not circularity: (a) Section 3 obtains the 1.25-1.5 yr intensity lag by interpolating a cross-correlation whose x-axis interval is 1 year to 0.01 yr; the CCF is defined only at integer-year lags, so the sub-annual peak location is not resolved by the data; (b) Section 2's trial-and-error top/bottom 5% window-intensity threshold is asserted to be insensitive to a few percent changes but no quantitative test is shown. Both could shift the claimed peaks/lags but do not make the analysis equivalent to its inputs. Because the central claim is a direct correlation of independent measurements with an external benchmark, the honest finding is no significant circularity (score 2 for minor, non-load-bearing self-citation).
Axiom & Free-Parameter Ledger
free parameters (3)
- intensity rejection threshold =
5% (top and bottom)
- contrast FWHM target =
0.10 to 0.11
- smoothing window =
5-point
axioms (4)
- standard math Cross-correlation is an appropriate measure of phase and amplitude relationship between annual time series
- domain assumption Supergranular network lane width is a proxy for magnetic flux at cell boundaries
- domain assumption Ca II K intensity at quiet Sun scales with network magnetic flux and temperature
- domain assumption The Kodaikanal 100-year data, after calibration, represent the true solar cycle variation without significant seasonal or instrumental drifts
Cite this review
Pith. "Pith review of Which Solar Latitude Follows the Sunspot Cycle Exactly?." pith.science (2026). https://pith.science/paper/UD6FZZSN
@misc{pith2026250909974,
author = {Pith},
title = {Pith review of: Which Solar Latitude Follows the Sunspot Cycle Exactly?},
year = {2026},
howpublished = {\url{https://pith.science/paper/UD6FZZSN}},
note = {Machine review of arXiv:2509.09974}
}
read the original abstract
The large-scale convection in the Sun known as supergranulation is manifested as a network structure on the solar surface. The network cells have an average lifetime of 24 hr, a size of about 30 Mm, and a lane width of about 6 Mm. We have obtained the lane widths and intensities at different latitudes from the Ca {\sc ii} K spectroheliograms from the 100 yr Kodaikanal archival data. We have then calculated the cross correlation function of lane widths and intensities with sunspot number at every latitude from 60$^{\circ}$ N to 60$^{\circ}$ S. The correlation coefficients of the quantities show an approximate North-South symmetry with broad peaks around $\pm$(11--22)$^{\circ}$ latitude with values of about 0.8. The results imply that these latitudes follow the sunspot cycle strongly. The maximum correlation for the lane widths occurs (18$\pm$2)$^{\circ}$ N and (20$\pm$2)$^{\circ}$ S with no phase difference. For intensities, this happens at (13$\pm$2)$^{\circ}$ N and (14$\pm$2)$^{\circ}$ S with a phase difference of 1.25 to 1.5 yr. It is interesting to note that the lane width correlations peak during the solar maximum whereas the intensitiy correlations peak 1.25--1.5 yr after the solar maximum. The results, generally show that no unique latitude exactly follows the solar cycle for all quantities. The results are important in flux transport on the solar surface and have implications for the quiet Sun UV irradiance variations.
Figures
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...
-
[3]
- [1] #1 = = ^ ^ ^ .\!\!^ d .\!\!^ h .\!\!^ m .\!\!^ s .\!\!^ @mss
thebibliography [1] 20pt to REFERENCES 6pt =0pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command E...
2021
-
[4]
Bappu , M. K. V. 1967, , 1, 151, 10.1007/BF00150312
-
[5]
Berrilli, F., Ermolli, I., Florio, A., et al.\ 1999, , 344, 965
1999
-
[6]
Chatterjee, S., Mandal, S., & Banerjee, D.\ 2017, , 841, 2, 70. doi:10.3847/1538-4357/aa709d
-
[7]
Chatzistergos, T., Krivova, N. A., & Yeo, K. L.\ 2023, Journal of Atmospheric and Solar-Terrestrial Physics, 252, 106150. doi:10.1016/j.jastp.2023.106150
arXiv 2023
-
[8]
Cossette, J.-F. & Rast, M. P.\ 2016, , 829, 1, L17. doi:10.3847/2041-8205/829/1/L17
-
[9]
De Rosa, M. L. & Toomre, J.\ 2004, , 616, 2, 1242. doi:10.1086/424920
-
[10]
Duvall, T. L.\ 1980, , 66, 2, 213. doi:10.1007/BF00150578
-
[11]
T., Phillips , K
Gallagher , P. T., Phillips , K. J. H., Harra-Murnion , L. K., & Keenan , F. P. 1998, , 335, 733
1998
-
[12]
Getling, A. V. & Kosovichev, A. G.\ 2022, , 937, 1, 41. doi:10.3847/1538-4357/ac8870
-
[13]
L., & Schou, J.\ 2003, , 421, 6918, 43
Gizon, L., Duvall, T. L., & Schou, J.\ 2003, , 421, 6918, 43. doi:10.1038/nature01287
-
[14]
Gizon, L., Birch, A. C., & Spruit, H. C.\ 2010, , 48, 289. doi:10.1146/annurev-astro-082708-101722
-
[15]
Gontikakis, C., Peter, H., & Dara, H. C.\ 2003, , 408, 743. doi:10.1051/0004-6361:20030940
-
[16]
doi:10.1007/s11214-023-01021-6
Hotta, H., Bekki, Y., Gizon, L., et al.\ 2023, , 219, 8, 77. doi:10.1007/s11214-023-01021-6
-
[17]
doi:10.1088/0004-637X/759/2/106
Huang, C., Yan, Y., Zhang, Y., et al.\ 2012, , 759, 2, 106. doi:10.1088/0004-637X/759/2/106
-
[18]
S., Bharati Das, S., Mani, P., et al.\ 2024, Nature Astronomy, 8, 1088
Hanson, C. S., Bharati Das, S., Mani, P., et al.\ 2024, Nature Astronomy, 8, 1088. doi:10.1038/s41550-024-02304-w
-
[19]
F.\ 1994, Solar Active Region Evolution: Comparing Models with Observations, 68, 1
Howard, R. F.\ 1994, Solar Active Region Evolution: Comparing Models with Observations, 68, 1
1994
-
[20]
Kariyappa, R. & Sivaraman, K. R.\ 1994, , 152, 1, 139. doi:10.1007/BF01473196
-
[21]
doi:10.1126/science.244.4901.197
Lean, J.\ 1989, Science, 244, 4901, 197. doi:10.1126/science.244.4901.197
-
[22]
Leighton , R. B., Noyes , R. W., & Simon , G. W. 1962, , 135, 474, 10.1086/147285
doi:10.1086/147285 1962
-
[23]
Lord, J. W., Cameron, R. H., Rast, M. P., et al.\ 2014, , 793, 1, 24. doi:10.1088/0004-637X/793/1/24
-
[24]
Mandal, S., Chatterjee, S., & Banerjee, D.\ 2017, , 844, 1, 24. doi:10.3847/1538-4357/aa76e3
-
[25]
McIntosh, S. W., Leamon, R. J., Hock, R. A., et al.\ 2011, , 730, 1, L3. doi:10.1088/2041-8205/730/1/L3
-
[26]
doi:10.1051/0004-6361:20030713
Meunier, N.\ 2003, , 405, 1107. doi:10.1051/0004-6361:20030713
-
[27]
doi:10.1051/0004-6361:20066644
Meunier, N., Roudier, T., & Tkaczuk, R.\ 2007, , 466, 3, 1123. doi:10.1051/0004-6361:20066644
-
[28]
doi:10.1051/0004-6361:20078835
Meunier, N., Roudier, T., & Rieutord, M.\ 2008, , 488, 3, 1109. doi:10.1051/0004-6361:20078835
-
[29]
H., Woehl, H., et al.\ 1989, , 213, 1-2, 431
Muenzer, H., Schroeter, E. H., Woehl, H., et al.\ 1989, , 213, 1-2, 431
1989
-
[30]
Patsourakos , S., Vial , J. C., Gabriel , A. H., & Bellamine , N. 1999, , 522, 540, 10.1086/307609
-
[31]
Priyal, M., Singh, J., Ravindra, B., et al.\ 2014, , 289, 137. doi:10.1007/s11207-013-0315-7
-
[32]
Priyal, M., Singh, J., Ravindra, B., et al.\ 2023, , 944, 218. doi:10.3847/1538-4357/acaf60
-
[33]
Raju, K. P.\ 2016, , 291, 3519. doi:10.1007/s11207-016-0991-1
-
[34]
P.\ 2018, , 478, 5056, 10.1093/mnras/sty1325
Raju, K. P.\ 2018, , 478, 5056, 10.1093/mnras/sty1325
-
[35]
Raju, K. P.\ 2020, , 899, L35. doi:10.3847/2041-8213/abacb7
-
[36]
Raju, K. P. & Singh, J.\ 2002, , 207, 1, 11. doi:10.1023/A:1015585010078
-
[37]
P., Singh, J., Ravindra, B., et al.\ 2023, , 959, L24
Raju, K. P., Singh, J., Ravindra, B., et al.\ 2023, , 959, L24. doi:10.3847/2041-8213/ad13e9
-
[38]
Reeves , E. M. 1976, , 46, 53, 10.1007/BF00157554
-
[39]
& Rieutord, M.\ 2018, Living Reviews in Solar Physics, 15, 1, 6
Rincon, F. & Rieutord, M.\ 2018, Living Reviews in Solar Physics, 15, 1, 6. doi:10.1007/s41116-018-0013-5
-
[40]
doi:10.1051/0004-6361/201423577
Roudier, T., S vanda, M., Rieutord, ., et al.\ 2014, , 567, A138. doi:10.1051/0004-6361/201423577
-
[41]
doi:10.1088/0004-637X/792/1/22
Sindhuja, G., Singh, J., & Ravindra, B.\ 2014, , 792, 22. doi:10.1088/0004-637X/792/1/22
-
[42]
Simon , G. W., & Leighton , R. B. 1964, , 140, 1120, 10.1086/148010
doi:10.1086/148010 1964
-
[43]
Singh, J. & Bappu, M. K. V.\ 1981, , 71, 1, 161. doi:10.1007/BF00153615
-
[44]
Singh, J., Priyal, M., & Ravindra, B.\ 2021, , 908, 210. doi:10.3847/1538-4357/abd021
-
[45]
Snodgrass, H. B. & Ulrich, R. K.\ 1990, , 351, 309. doi:10.1086/168467
doi:10.1086/168467 1990
-
[46]
Sridharan, R , 2017--18 IIA Annual Report
2017
-
[47]
S \'y kora, J.\ 1970, , 13, 2, 292. doi:10.1007/BF00153550
-
[48]
Tian , H., Marsch , E., Tu , C. Y., Xia , L. D., & He , J. S. 2008, , 482, 267, 10.1051/0004-6361:20079235
-
[49]
Wang, H.\ 1988, , 117, 2, 343. doi:10.1007/BF00147252
This paper was first reviewed by deepseek-v4-flash on August 4, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.